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ABSTRACT 



Under the assumption of constant failure rates it is pos- 
sible to build a "reliability shorthand" which gives a sim- 
ple, unified approach to reliability computations for sys- 
tems in the presence of complications like support by shared 
spares or changes in the failure rates of surviving compo- 
nents when other components fail. The computational imple- 
mentation of the shorthand depends upon the convolution of 
strings of exponentially distributed random variables. 

This paper presents an algorithm for the numerical convo- 
lution of exponentially distributed random variables. After 
reducing the system scenario to its shorthand format, one 
can use the programs that are given in the appendix to 
obtain numerical values for the reliability of the system. 
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I. INTRODU CTJON 



The reliability shorthand considered in this paper has 
beer, developed in conjunction with the course OA 4302 "Rsli 
ability And Weapon Systems Ef fectiveaess Measurement” at th 
Naval Postgraduate School. A tutorial introduction to the 
reliability shorthand was given by Repickyf Ref . 2 ]. This 
paper is devoted to a complementary part of the idea. 

Any study on system reliability always requires -.wo 
steps; one is the description of the system’s life and the 
other is the derivation of its survival function. 

Under the assumption of constant component failure rates 

♦ 

this paper presents a way of obtaining the system’s reli- 
ability which requires little beyond the description of the 
system's life. 

Section II contains an approach to the convolution of 
exponentially distributed random variables. Also there is 
presentation of a computational algorithm for the convolu- 
tion of exponentially distributed random variables. 



7 



Appendix A gives survival functions corresponding to sev- 
eral reliability shorthand norations and a program writer in 
Fortran for computations from shorthand notations. 

Section III deals with the reliability of redundant sys- 
tems under the assumption of constant component failure 
rates . 

Appendix B consists of a program written in Fortran. The 
program supports the approach of Section III. There is a 
crude Monte Carlo simulation program in Appendix D which is 
a simulation program parallel to the program in Appendix 3. 

There is another program in Appendix C written in For- 
tran. This program uses the network approach to systems 
described in Section III. 

Appendix E summarizes the definitions used in this paper. 
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II 



AN APPROACH TO COMPUTING CONVOLUTIONS OF EXPONENTIAL 

m ggrmiAiiis 



This section introduces a general algorithm for com- 
puting the survival function of any convolution of exponen- 
tial random variables. 

In reliability, the term convolution refers to rhe summa- 
tion of independent random lives. In order to have simplic- 
ity in specifying convolutions, the reliability shorthand 
introduces a special notation. In the following sections we 
will use this notation. 

A. THE SURVIVAL FUNCTION FOR A CONVOLUTION OF RANDOM 

VARIABLES 

Let F 1 (t ) and F2 (t) be the survival functions for the 
random variables T1 and T2 respectively. Let fl(t) and f 2 (t) 
be the corresponding densities. Let "F (t) be the survival 
function for the random variable T, where T=T1+T2. 

Then the likelihood expression for "F (t) is 
t 

F (t) =F1 (z) + / F2 (t-s) fl(S) ds 
0 



9 



In the right hand side of equation, FI (- 1 -) is rhe prob- 
ability that component one completes the mission, 
t 

F 2 ( t-s) f 1 (s) ds is the probability that at some time 
s (0<s<t) component one fails, component two takes its place 
and carries rhe system to the end of the mission duration t. 

In order to illustrate consider soma applications. 

1 . Example 

Reliability Shorthand Notation : EXP {Xi } +EXP (X 2 } 
SYSTEM: One component having one spare with a dissimilar 

failure rate. If the active component fails, the 
spare will replace it immediately. 

Here the life for the system is T = T1 + T2. The reli- 
ability shorthand notation indicates that this system has an 
exponential life with failure rate Ai followed by an expo- 
nential life with failure rate A 2 . 

The survival function for the active component is 

f"-, (t)= e‘ Xlt , t>0. 

The survival function for the spare is 
? 2 (t) = e‘ X2t , t>0. 

The survival function for the system is 
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F(t) = E (t) ♦ / F ? (t-s) f n (s) 5s 

F(t) = e ‘ Xlt + J e " x 2^ t_s ) Ai? XlS ds , t>0. 

If we complete the integration, the result: is 



F (t) = \ 2 /(A 2 - Xi) e Xlt + Ai/ (X!-X 2 ) e" X2t , t>0. 



Which is the well known result. 



Another way to establish this formula is the use of the 
moment generating function. (Freund and Walpole[ Ref . 3 ]) 

2 . Ex ample 

Reliability Shorthand Notation : EXP { A } + EXP { A } 
SYSTEM: One component having one identical spare. If the 

active component fails, the spare will replace it 
immediately. 

The formula that we derived above becomes meaningless, 
because the denominators become zero. If we proceed as 
before 

t 

F(i) = F1(t) + J* F2 (t-s ) fl (s) ds 

0 




0 



F (t ) = ( 1 ♦ A*- ) 5 At , t>3 . 
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The result comes out as expected co be the Erlang (2, A } 
survival fur.cxion. 



B. RELIABILITY SHORTHAND NOTATION 

EXP {A} + EX? {A} +. . . . +EXP{ A} (n identical 

exponential lives) 

This leads to the Erlang {n, A } survival function 

_ n , . 

F(-) = l (At) /(i-1) ! e~ At , t >3 . 
i =1 

The use of moment generating function gives the result 
immediately. 

C. THE RELIABILITY SHORTHAND NOTATION 

EXP (Ai) +EXP {' A 2 } + + EXP (A n ) 

This is the expression for the convolution of n random 
variables where each has a distinct failure rats. 

By the approach used in Section 2.1, adding one exponen- 
tial life at once, one can obtain the formula for the survi- 
val function 
n 

F (t ) = l n. A / n ( A j ”A 2 ) e’V , t>0. 

1=1 j*l J j.*l 
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This is also a w^ll known formula and can be obtained 



from by a moment generating function. 

D. THE RELIABILITY SHORTHAND NOTATION 

EXP {Ai } + EXP (A 2 } + + EXP{A 2 } (n t terms) 

+ EX? {A 2 }+EXP {A 2 } + +EXP { A 2 ) (n 2 terms) 



+ EXP[A k } + EXP{A k } + +EXP{A k ) (n k terms) 

k 

This is the convolution of F n. exponential random 

i-i 1 

variables where there are n^ identical exponential random 
variables having the failure rate A.. . 

The moment generating function technique is not useful 
this situation, since there is a huge amount of complexity 
involved. This section deals with this notation using the 
convolution formula. 

1. Rel ia bility Shorthand Notation 

EXP { Ai > +EXP { Ai } + EX? (A 2 ) 

We know the survival function for EXP {Ai } +EX? { A 2 } and 
also for EX? { A }+EXP{A). The next two subsections present 
different ways no reach the survival function for the shcr 
hand notation above. 
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(a) 



Let T1,T2,T3 be random variables distributed as 
EXP {Xi } , EX?{Ai}, EXP{A 2 ) respectively. 

Let T1'=T1+T2. This random variable has the Erlang {2, A} 
distribution. 



The convolution formula for T=T1'+T3 is; 



t 

_ 



F T (t)=F T/ (t) + 
F r (t) = (1 + A : t)e' Xlt + 



F Ti (t-s) f T (s) ds , t>0 . 
t 1 



e- XlS ds 



-A 2 t 



F (t) — e [ ( A 2 2 -2 XiX 2 ) / (A 2 -A 1 ) 2 + AxA 2 /(A 2 ~A i ) 

T 

+ e" X2t { A 1 2 /(x 1 -X 2 > 2 } ,t>0. 



(b) . 



Let T1,T2,T3 be random variables distributed as 
EXP {Ai} , EXP {Ax}, EXP{A 2 } respectively. 

Let T2’=T2+T3. This random variable has the survival 



function 



- Aj t 



r A 2 t 



F (t) =A 2 /(A 2 -Ax)e + Ax/(Ax-A 2 ) e A2 

‘2 



, t >0 . 



The convolution formula for T=T1+r2' is 



F j (t) = F t / (t) ♦ 



F t (t-s) f T ,(s) ds , t>0 . 
1 1 < 2 
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f t (t) 



-Ait 



-x 2 t 



— X2/(X2~Xi) S + Xi/(Ai~A.2} 9 + 

( e“ XlS -e ^ s ) ds 



iXl liA^(A 2 -A 1 ) 



— -A 1 t 

Fj (t) = e C(A 2 2 -2 X 1 X 2 ) /(A 2 -Ai > 2 + Ai /(A 2 -A i) t} 

*~e zt {\f / ( A x “ A 2 ) 2 } » - — 0 • 



Subsections (a) and (b) illustrate that the convolution 
formula gives a unique result, regardless cf the way of 
choosing the prior random life. 

2. Re liabi lity Shorthand Notation 
EXP {Xj} + EXP (Ai } + EXP {A 2 } + SXP {x 2 } 

Let T 1, T2, T3, T4 be random variables distributed 
as 

EXP{Ai}, EXP { Ai } , EXP{A 2 }, EXP{x 2 } respectively. 



Then from the derivation in Section 2, T1 '=T1+T2+T3 is a 
random variable having the survival function 
F t / (t) = (a +a t) "e Xlt ta 2 i€ Al H^0. 

I.j* 11 12 

where an = (A 2 2 “ 2 AiA 2 ) / (X 2 “ Ai ) 2 • (X 2 “Xi ) • 

a2i = Ai 2 / <Ai“A 2 ) 2 • 



The convolution formula for T=T1'+T4 is 

- f- 

F t (t) = F-^t) + J F (t-s) ^/(s) ds , t>0. 



_ -Ait -A 2 t 

F T (t ) =(an+ai 2 t)e +a 2 ie + 

+ A 2 a 2 je 2 ) 3 s 



e Az ( t jAj-aj, 2 +A 1 ai 2 s) e lS 
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t 




The result from the above is 



F T (t)= (a 1 ‘ 1 +a 1 '^)e Xlt +(a 2 ^^ 2 t> e' Alt ,t>0. 
where a x V= (Xz>-3 A^Aj ) /(A 2 - Xj) 3 , V 2 * ^iA 2 2 / (A 2 - Xi ) 2 / 

a 2 \= Ui 3 -3A 1 2x 2 )/(X 1 -A 2 )3 e h , 2 =X 1 2 X 2 /(X 1 ~X 2 ) 2 • 

It is important to none that the number of exponential 
nerms equals the number of dissimilar failure rates and each 
exponential term has a polynomial coefficient with the 
degree of the polynomial equal to the number of identical 
random variables having the corresponding failure rate. 

The next section deals with the convolution of exponen- 
tially distributed random variables using the fact illus- 
trated above. 

3. In trodu ction of an Algorithm for the Convolution of 
Exponential Lives 

From Subsection 2, we can infer the form of the sur- 
vival function which we sought at the beginning of 
Section C. 

SHORTHAND NOTATION : 



EXP { X i} + EXP [Xi } + . .. 




(t L -.eras) 


EXP {A 2 } + EXP (A 2 } + ... 




(n t terms) 


EXP(A k }+2X?{A k }+... 


... +EXP(A k J 


(n^ o°ras) 
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SURVIVAL FUNCTION : 



F(t) = li(t)e^ 2t + A 2 (t) 9 A2t +...+ A (t) e A k t , t >0 

n, -i 

where Ai (t ) =a 1 1 +a 12 t+ a t 2 *....+a 2n t , 

1 3 1 

n 2 -i 

A 2 (t) =a 21 +a 22t + a 23t 2 + . . . . +a 2n t , 

1 



A k ( t )=a kl +a k 2 t+a t 3 t 2 + 



+ a 



kn. 



tV 1 



a. Example: 

SHORTHAND NOTATION : 

EXP { A'i } + EXP {A 2 } 

+ expca 2 } 

+ EXP{A 3 }+EXP{A 3 }*EXP{A 3 } 

SURVIVAL FUNCTION : 

F (t) = { a 3 1 +a 12 t) e 1 +a 2 3 e 2 + {a 3 i+a 32 'C+a 33 t 2 )e^, "t>0. 
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b. An Algorithm to Compute the Coefficients 

The algorithm represented below develops the 
survival function by adding one random variable in each run 
As an example , in order to compute the survival function 
for the convolution of ten exponentially distributed random 
variables, the algorithm is supposed to run ten times. 

The notation used in the algorithm is: 



K 


number of dissimilar failure rates 




ith type failure rate 


A. 


failure rate for the currently entering life 


a Jk 


kth coefficient on the jth polynomial 


n i 


current number of identical lives having 
the ith failure rate 


nn. 


number of random variables having i*h 
failure rate. 


Initial: 


aj k =0 , v where j=1,2,..K , k = 1,2,...r.j 

n.j=0 , v- where i= 1 , 2 , . . K 


Input : 


A^ , v- where i=1,2,...K 
nn.j , V i where i = 1 f 2,...K 


The firs 


t run is : n -j= 1 , a-j -j = 1 
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Algorithm : 



n.=n.-+1, until n.- =n n • . 

'e 'e ’e 'e 

1. Update the coefficients : a. ^ for k=2,3,..,n. 

e e 

a k = A . a • / (m-1) , where m=n. -j for j=0,1,.,n. -1 

'e K e e m-1 ’e 'e 

2. Update the coefficient : a.. ^ 



n . 



3i =a- !+ I (a i 1 Ai/(A^ i )+ l ( j-1 ) ! A. * . / ( a 2 -A.) J } 
e e i?H e 1 e j=2 e J \ 



^.>0 



V 



3. Update the other coefficients: a.^ Vi»k 

where i*i , n.*0 for i=1,2,.,K, k=1,2,.,n. 

e l i 



ai "t' a i n iV < V Xl1 ’ Vi " h ' Ce ^ ' 

a in =(A i a in *" s imH )/<X f'V * vi “ h ' re iili e aad r -i-’ 
e e 

for m=n^-j, j=1 , 2 , . . , n.. - 1 



c. Example: 

Reliability Shorthand Notation : 
EXP&i } + EXP( Al }*EXP{A 2 }*SX?a 2 } 



Let T1,T2,T3,T4 be random variables distributed as 
EXP {Xi } r 2X? (A ! } f EXP[A 2 }, EXP{A 2 } respectively. 
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We would like to derive a formula for the survival func- 
tion of the random variable T=T1 +T2+T3+T4. 

The use of algorithm: 

1 st RON; 

n x = If i = 1 

At the end of this run we have only one random variable, 
which is distributed EXP{A 1 }. 

2 nd RON; 

i = 1 , n, =2 
e - 1 

a i2 ~ Ai&i i then s 12 = Ai 

a ii =a ii + 0 then a n = 1 

At the end of the 2. run we have the survival function 
for the T’=T1+T2, where T1 and T2 are identically distri- 
buted as EXP £Xi } • The survival function is 

F(t) = (a 11 +a i2 t)e‘ A ^ 

F (t) = ( 1+ Ait) e” Al , t t> 0 

Which is ERLANG {2, Ai} . 

3 ^ RON; 

i =2, n =2, n =1 
e 1 z 
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a 2 i =0+ai 1 A 1 /(A 1 -X 2 )+a 1! A 2 /(Ai-A 2 ) 2 



*1Ai/ (Ai- h) + AiA 2 / (Aj -A 2 ) 2 
= Ai 2 / (Ai“A 2 ) 2 

a i2 =a 12 A 2 /(A 2 -Ai) then a 12 = AiA 2 / - x 1 ) 

an =(A 2 a ir aj 2 ) / (a 2 “Aj then a ia = (A 2 3 ~2 AiA 2 ) / ( A 2 -A a ) 2 
Note that, here the coefficient a 12 is the updated one. 

At the end of the 3. run we have the survival function 
for the random variable T”= T1+T2+T3 , where T1,T2,T3 are as 
defined before. 

“Ait -A 2 t 

F-jiift) = (aij*-a 12 t) e +a 22 e , t>0 

where %,= (a 2 3 -2x 1 A 2 ) / (A, "Ai ) 2 . **= AiA/<A 2 "A 2 ) ' 
a 2 r= Ai (Ai ”A 2 ) 2 • 

Note that the coefficients are identical to the result 0 
Subsection (b). 

4 th RUN; 

i =2, r-j =2, n =2 

e z 

a 22 ~ A 2 a 2 j th en ~ 22 ~ ^1 ^ A 2 /(Ai”A 2 ) 2 

a 2i — a 2i + ^iiAi/(Ai ”A 2 ) 5 ^ 2 / (a i-A 2 > 2 

then a 21 =(x 3 *3 a 2 2 A 2 ) / <Ai “ A 2 > 3 
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a i 2 =a i 2 A 2 /(A 2 -Ai) then a 12 =AiX 2 2 /(A 2 -\i ) 2 

M = ^A 2 a n -a 12 )/ (x 2 -Xj) then = (X 2 3 -3X 1 X 2 2 ) / (X 2 '\ 1 ) 3 

At the end of the 4. run, we have the survival function 
for the variable T=T1 +T2+T3 +T4. 

F T (t) = (a!#a 12 t) e Alt +(a 2I +a 2 ^) e X2t ,t>0 

where an =(X 2 3-3A 1 A 2 2 )/(A 2 'Ai) 3 * a i 2 = A x A 2 2 / (a 2 * a x > 2 ' 
a 2 i = (A 1 3 -3A 1 2 A 2 )/(A 1 -A 2 ) 3 ' a 22 = Ai 2 A 2 / ( Ai ~A 2 ) 2 • 

Note that the coefficients are identical to the result of 
Subsection 2. 
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III. RELIABILITY shorthand approach TO SYSTEM RELIABILITY 

This section deals with a system whose components have 
constant failure rates. 

Having the reliability network for a system in which each 
component has an exponential life and knowing the probabili- 
ties for failures (discussed in Appendix E2) makes it easy 
to describe the system's life. 

In order to make the idea clear, we will go through some 
examples. 

A. EXAMPLE 

A 1 A 2 




2 spares ( Ai ) 1 spare (X 3 ) 

SYSTEM : 2 components in series with spares. 

DESCRIPTION : System has 2 components in series. Compo- 
nent 1 has a life distributed as SXP{Ai} and there are 2 
identical spares. Component 2 has a life distributed as 
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15 XP { A 2 } and there is one non-identical spare, whose life is 
distributed as EX?{A 3 }. 

LIFE : 



yEXP{Ai+A 2 } 

3 Ai/(A! + A 2 ) 

^EXP£ Ai + A 2 }+^ 

^ Aj/(Aj + A 2 ) ^ A 2 / (Ai + A 2 ) 

/ \ 

EXP{A 1 + A 2 }+^ EXPfAj+A,}* 

2 X 2 / (A 1 + A 2 ) 3 a 3 /(Ax + A 3 ) 

\ 

/ZERO ZERO 



2XP{A! + A 3 } 

7 Ai / (Ai + A 3 ) 



^ 1 ,/ h , +) 



6 A 3 /(Ai + A 3 ) 



EXP JAi+A^ ♦' 



5 Ai/(Ai+A 3 ) 

\ / 

x Ex?( Al n 3 }♦, 



yEXP {Ai + A 3 } 
9 Ai/ (Ai + A 3 ) 



\ 

10 A3/(Aj + A 3 ) 

\. 



‘ZERO 



Explanation for the derivation of system’s life: 

At the beginning the system has an EXP{Ai+A 2 } life. The 
failure of the system is by the failure of component 1 with 
a probability of Ai/(Ai+A 2 ) or by the failure of component 2 
with probability of A 2 /(Ai +A 2 ). 
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In the life figure, number 1 denotes the event "failure 
of component 1" and number 2 denotes the even 4 : "failure of 
component 2". 

If event 1 occurs, component 1 is replaced by one of the 
spares and system again functions with a life distributed as 
EXPtXj+Xzl# since the exponential distribution has the memo- 
ryless property and component 2 still has the same failure 
rate. 



If the event 2 occurs, 
spare and the system has a 



component 2 is replaced by its 
life distributed as EXP(Xi + X 3 } . 



The numbers on the life figure correspond to the transi- 
tions that can occur. The probability on each arc shews the 
conditional probability of the transition. As an example 
event 3 can occur with probability of Xi /(Xi + X 2 )* given that 
event 1 has occured before. 



The distribution ZERO defined by ssary[Ref.1] and 
Repicky f Ref • 2 ] (also defined in Appendix Z3) enters when 
life is exhausted. 

For convenience of description, it is helpful to define 
the concept of path used in this paper. Path denotes the 
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sequence of events in the system's life from the starring 
point to the point where the system is not functioning. 

Exa mple s ; 

Events 1 and 3 are a path, which denotes a sequence of 
lives for the system. In this case, the system has 3 expo- 
nentially distributed lives EXP {A a + A 2 } +EXP {\ 1 * Xl * +5XP C\i + A 2 * 
with the probability of Ax /( Xl + A 2 > A 1 /(a 1 + A 2 i * 

Events 1, 4 and 8 form another path, which describes a 

sequence of lives for the system. In this path, the system 
life is EXP { A ! +A 2 } +2 XP(Ai + A 2 } +EXp fAi+ A 3 } +ZERO. The prob- 
ability of this path is 

[Ai/(Ai +A 2 }(a 2 / (A i ♦ A 2 ][a 3 / (Aj +A 3 jj • 

The ZERO distribution contributes zero additional life to 
the system, so we can omit it. Nevertheless, we can not omit 
its probability in the calculation of path probability, 
since the event numbered 9 has a probability of occurring. 
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B. THE SURVIVAL FUNCTION FOR THE SYSTEM 



In this section we will deal with an approach to obtain- 
ing the system’s survival function. 

1 • EXAMPLES 

SYSTEM: Two components in parallel. 



Ai 




LIFE : 



M /Ui + X*) SXP{a 2 } 



EXP{Ai + A2} + 



X 2 /(Ai + x 2 ) 2 XP { A x } 



DESCRIPTION: Component 1 has EXP (a x } life and component 2 
has EXP(A 2 } life. 



There are two failure events that can be allowed. The 
first is that component 1 fails at some time t and component 
2 carries the system for the rest of the time. The ether one 
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is that component 2 fails at some time during the mission 
and component 1 carries the system for the rest of the time. 
There are more events such as no failures during the mission 
duration, which are taken care of by the representation. 

Now we have two paths 

No . of Path We ight Life 

1 Pi = Ai/(Ai + A 2 > SXPUj+AzJ+EXP^ } 

2 p 2 =a 2 /(a 1 + a 2 ) exp{ Ai + A2 } + exp{ Xi } 



Let T be the system's time to failure and let T1,T2,T3 be 
random variables exponentially distributed with the failure 
rates Ai+A^Aj./Aj respectively. 

Then 



' T 1 + T2 



T = < 



T 1 + T3 



with probability Pl 



with probability p^ 



The survival function can be written as 
F(t) = Pl F l (t) + p 2 ? 2 (t) , t>0 

where Pl = Ai / (\i+\ 2 ) , P 2 = A 2 /(Ai + A 2 ) and F x (t) , F 2 (t) denote 
the survival functions for the shorthand notations 
EXP {A x +A 2 } +EXP [A 2 } » SXP(Ai ♦ A 2 } *2X? {\ x } respectively. 
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The survival function for the convolution of two exponen- 
tially distributed random variables with dissimilar failure 
rates is 

F(t) = A 2 / (A 2 - Ai) i Xlt +Xi/(Ai-A2 ) 9 A2 t>0. 

If we do the substitutions for Fi (t) and F 2 (t) as 
A i =Ax +A 2 , A 2 = A 2 andAi = Ai+ A 2 , A 2 = Ax respectively, 

Fj (t) becomes 

'? 1 (t) = A 2 / (- Ax) e^ Al+A2 ^ + (Ax +A 2 ) /Ax e X2 ° , t>0, 
and F 2 (t) becomes 

-(Ax+A 2 )t -A 2 t 

F 2 (t) =A x /(-A 2 ) s + (Ax + A 2 ) /A 2 £ / t>0. 

Then 

F(t) = Ax/(Ax + A 2 ) F x (t) + A 2 / (Ax + A 2 ) ? (t) 

F(t) = e Al +e Al 4 Al+Az)t , t >0 . 

The result gives the survival function that is well known 
for this system. 

Another example is the 2 out of 3 system. 
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RELIABILITY NETWORK: 




The known survival function is 

■ ( A. i +X 2 ) t -(A 2 +A 3 )t - ( A^+A 3 )t -( Ai+A 2 +A 3 ) t 



F(t)=e 



+e 



+ s 



-2s 



, t>0. 



LIFE: 



,EXP{A 2 +A 3 » 



Ai / ( Ai + A 2 A 3 ) 



EXP {Ai +A 2 + Aa } +— A 2 / (Ai + A 2 + Aa ) — EX?{a x + A 3 ) 



A3/ (Ai + A 2 + A 3 ) 



'EX P (Ai + A 2 } 



ZX? [ Ai + A 2 +A 2 } +EXP {A 2 + A 3 } with probability \y (\f\ z +\ 3 ) 



T= \ EXP { A x + A 2 +A 3 } +EXP {Ai + A 3 } with pro babilit y a 2 / ( 1 !%%) 



EXP f Ai + A 2 + A3 } +EXP {Ai+A 2 } With probability \/ (a 
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The survival function is 



?(t)=Pj Fx (t) + p 2 F 2 Ct)+P3 F 3 (t) , t>0. 



where 



-(Ax+A 2 +X 3 )t ~(A 2 +A 3 )t 

(t) = (A 2 +Ai/ (”Ai) e + (Ax<A 2 +X 3 ) /Ax e 



?z (*) = (Ax+A^ /(-A 2 ) e 



- C A x + A 2 + A 3 ] it 



(AaVxP Aj 



■CAx+A 3 ) t 



— -(Ax+A 2 +X 3 )t -( Ax+A 2 ) t 

F 3 (t) = (Xx+X>) / (*X 3 ) e + (Ax+Xz+Xa) A 3 e 

If we do the necessary cancellations, we can get, 

p (t)=e - (A 2 +A 3 )t +e -(Ax+X 3 )t +e -(Ax+A 2 )t _ 2e -(Ax+A 2 +A 3 )t } t>Q 



as desired. 



2. General Procedure 



Having the algorithm presented in Section II, we 
treat mere complicated systems under similar assumptions 

The procedure is 

i. Set up the reliability network for the system. 

ii. According to this network, set up the system life 

iii. Using the proper reliability shorthand formulas 
the related convolutions of exponential lives, set up 
survival function. 



can 



for 

the 
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IV 



StmAHI 



The reliability shorthand is an easy way to describe a 
system’s life, but it is difficult to implement computet io- 
naly since there is considerable complexity in handling con- 
volut ions. 

The algorithm presented in this paper gives some relief 
from this difficulty. However, the accuracy in obtained from 
this algorithm is very much related to the differences in 
the failure rates. 

Another aspect in the algorithm is that distributions are 
convolved one at a time and this requires very accurate run- 
ning conditions in the case of a complicated system. 

It is believed that it is possible to derive another 
algorithm which is more powerful than the one introduced 
here. Instead of adding one distribution at a time, one can 
try to add several at a time. 
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APPENDIX A 



A- 1 

This section contains the survival functions for several 
shorthand notations which were derived by the use of the 
approach described in Section II. 

A. 1.1 Shorthand Notation : SXP{A}. 

Survival Function: F(t)= e At , t>0. 

A. 1.2.1 Shorthand Notation: SX?{A!}+ SXP{A 2 } 

Survival Function; 

F(t) = A 2 /(A 2 -Ai) e Xl ♦A 1 /(A 1 -A 2 ) e A2t ,t>3. 

A. 1.2.2 Shorthand Notation : EXP ( A } + EXP { A } . 

Survival Function: 

F (t)= (1+ A t) e' At , t>0. 

A. 1.3.1 Shorthand Nonation : EXP } + EXP {x 2 } + ZX? } 
Survival Function: 

F ( + ) =an e" Alt + a 21 e A2t *a 31 e Aat , t>0. 
where a x x = A 2 A 3 / (a 2 * Ax ) (a 3 ”A 1 > ' a 21 = XlA3/ (Ai -A 2 ) (A 3 -A 2 ) , 
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a 3 1 — X i A 2/ ( A 1 X 3 ) ( X2 "" X 3 ) - 



A. 1.3.2 EXP{A } + EXP{ A } ♦ EXP{ A } 

Survival Function: 

F(t) = (1+A t+1/2A2t2) 9 ' Xt , t>0. 

A. 1.3.3 EXP{Ai } + EXP{A 2 } ♦ EXP(A 2 } 

Survival Function: 

F(t) = au e’* lt + (a +a t)e"^ zt , t>0 

21 22 

where a n * A 2 2 / ( A 2 -A x ) 2 , hi~ tti 2 -2 a x A 2 ) / (A x - A 2 ) 2 
^2 2f~ \\\?/ ( A 1 *A 2 ) • 

A. 1.4.1 EXP {A 1 } +EXP{A 2 } <■ EXP ( A 3 } +EX? {A 4 } 
Survival Function: 

F(t) =aue Al +a 2 ie A2t +a 3 re X3t +a f re X4t , t>0. 

where a.-, = TI Aj/ n (A. -A.) V i=1,2,3,4 

1 ' J>i J>1 J 

A. 1. 4.2 EXP {A } +EXP [ A } ♦ EXP { A } + EX P [ A } 
Survival Function: 

F (t) = ( 1+At+1/2 A 2 t 2 +1/6 A 3 t3)i At , t>0 . 

A. 1. 4.3 EXP {Ai } +EXP(A 2 } + EXP(A 2 } + EXP {A 2 } 
Survival Function: 
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t>0. 



-Xit -x 2 t 

F(t)=a il « + ( a 21 +a 22 t + a 23 t 2 ) e ' 

where a x = X 2 3 / ( X 2 -X x ) 3 , \ 1 =1 ~ a n 

a 22~ X x X 2 (X “2 X 2 ) / (a X 2 ) 2 • a 23 = X x X 2 ^A x ~X 2 ^ * 

A. 1.4.4 EXP[\i }+EXP(Xi} *EXP{A 2 }+EXP[A 2 } 

Survival Function: 

— -X2t -X 2 t 

F (t ) = ( a il +a i2 t ) e + ^ a 21 + ^2 t ) 5 ' t -°* 

where a : f ( A 2 3 -3a 2 2 Xi ) /(X 2 " Xi ) 3 , a 12 = X x X 2 2 / (X 2 - A x ) 2 r 
a 2 i = ( X 1 3 -3xfX 2 ) /(X x -X 2 ) 3 , a zrXzXiVUi-Xz ) 2 . 

A. 1.4.5 2X?(Xi } +EX?(X 2 } *• EXP{X 3 } + EX? {X 3 } 

Survival Function: 

-Xit -X 2 t -Xst 

F(t)=a e +a e +(a + a t) e , t>0. 

11 21 3 1 32 

where a li = ^2X 3 ^(x 2 -a x ) (X 3 -X x ) z , a a = X x A 2 2 /(A x “A 2 ) (X 3 -Xj> ) z • 

a 3 ~ Xi X 2 / ( Xi ” X 3 ) (X 2 -X 3 ) + X x X 2 X 3 [ 1 / ( Xi “X 2 ) (A x ~ A 3 ) 2 

~V(Xi *x 2 ) (x 2 - X 3 ) 2 3 

a 32 = ^ 2X3 /( X a ” X 3 ) (X 2 -X 3 )* 
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A- 2 



This section introduces a Fortran program using the 
algorithm described in Section II. 

A PROGRAM FOR THE ALGORITHM ONE AT A TIME: 



THIS IS A PROGRAM TO COMPUTE THE RELIABILITY OF A SYSTEM 
WHICH HAS THE RELIABILITY SHORTHAND NOTATION 
EXP (LI) . +EXP (LN) 

WHERE THERE IS NO RESTRICTION FOR THE FAILURE RATES. 
VARIABLES: 

A (If J) : I.TH TYPE FAILURE RATE, J.TH COEFFICIENT 

ON THE POLINOM. 

NI (I) : AMOUNT OF LIVES HAVING THE I.TH TYPE 

FAILURE RATE. 

NINIT(I) : AUXILAURY ARRAY FOR THE NI (I) . INPUT 

FOR THE PROGRAM. 

L (I) : THE ARRAY FOR THE FAILURE RATES 

REAL L(20),A (20.10) 

INTEGER NINIT (20) , NI ( 20) 

GET INPUT 



2 

4 

6 

9 

3 



READ (5,49 9) T 
WRITE (6 ,498) T . 

READ (o , 50 0) K 
WRITE (6, 1500) K 

READ (5,501) (NINIT (I) ,L (I) ,1=1, K) 

WRITS (6,1501) ( NINIT ( I) ,L (I) ,1=1, K) 

COMPUTE COEFFICIENTS ONE ALL 
DISSIMILAR 



KK=1 



A (KX, 1) =1 . 
NI (KK) =1 

JJ=1 



IF { JJ . EQ . KK) 30 TO 6 
A (KK, 1) =A (KK, 1) *L (JJ) /(L(JJ) -L (KK) ) 
IF(JJ.EQ.K) GO TO 9 
JJ=JJ+1 



GO TO 4 

IF (KK.EQ. K) GO TO 8 

KK=KK+1 

GO TO 2 

CONTINUE 



BEGIN TO ADD ONE AT A TIME 



I E=1 

32 IF (NINIT(IE) . EC .NI (IE) ) GO TO 99 
NI (IE) =Nr (IE) +1 
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UPDATE IE 



nnn non nnn 



12 



20 SUM=0. 



J = 0 

N NNN=NI (IE) - J 
A (IE, NUNN) =L (IE) *A 
IF (NNNN.EQ.2) 30 I 
J = J+1 
GO TO 12 



(IE, NNNN-1) /FLOAT (NNNN-1) 
0 20 



1 02 

101 



21 



DO 101 1=1, K 
IF(I.EQ.IE) GO TO 101 
SUM=SUM+A (I, 1) * 

IF (NI (I) . LT. 2) 

FACT=1. 

NKK = NI (I) 

DO 102 II 
FACT=FACT 
SUM=SUM+ A (I ,11) * F 
CONTINUE 
CONTINUE 
A (IE, 1) =S UM+ A (I 



G^TO - L(IE,) 



= 2 ,NKK 
<II-1) 

ACT* L (IS) / (L (I) -L (IE) ) **II 



E, 1) 



UPDATE I. NS. IE 



24 

26 

99 



1=1 

IF(I.EQ.IS) GO TO 26 
A(i,Nim )=A(I,NI(I> ) 
IF (NI jl) . EQ. 1) GO 

NNI = NI fl) -J 
A (I , NNI) = (L 
* / 

IF (J. EQ. (NI 
J=J + 1 
GO TO 24 
EQ. K) GO TO 32 




*L(IE) / ( L ( IE) -L (I) ) 
TO 26 



*k (I, NNI) -FLOAT (NNI) *A (I, NNI+1) ) 
E) -L (iy) 

1)) GO TO 26 



IF (I 
1 = 1+1 
GO TO 21 
CONTINUE 
IF (IE . EQ. K) 
IE=IE+1 
GO TO 32 



GO TO 35 



CALCULATE THE PROBABILITY 
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41 

40 



CONTINUE 

PRO=0. 

DO 40 1=1 , K 
S U M= 0 . 

NKK= NI ( I) 

DO 4 1 J=1 , NKK 
SUN=SUM+ A (I, J) *T*< 
CONTINUE 

PRO= PRO + SUM*EXP (-L (I) 
CONTINUE 



J-1) 

*T) 



PRINT OUT THE RESULTS 



17 



498 

499 

500 

501 



WRITE (6,600) 
DO 17 1=1 ,K 
NNK= NI ( I) 
WRITE (6,1600) 
WRITE (6,603) 
WRITE (6,601) 
STOP 
FORMAT 
FORMAT 
FORMAT 
FORMAT 



L (I) , NI 



ii’i. 



NNK) 



M«,10X,'THE INPUT IS : » , // 1 5X , » TI ME= » , ? 1 0 . 4) 
FI 0 . 4) 

3X,I3) 

3X , 13, FI 0. 5) 



37 



600 POEM AT ( 1 1 ’ , 1 OX, ’FAILURE RATS 
*' ,/1 OX, 40 (*-• ■ 

6 01 POSH AT 



*F10. 7) 



'('OMol 1 , 'PROBA BILITY ( TIME > ’,F10.5,') = 



AHOONT COEFFICIENTS 



6 03 POEM AT ( '0‘ ,35X, 3 (5 (El 2.5,2X) ,/) ) 

1500 FORMATj'O' , 14X, »K = «, 3X , 12 , 7X , » (THE NUMBER OF 
♦DISSIMILAR FAILURE RATES) *) 



1501 FORMAT (’O’ .15X, '# OF R.V. 

125 ( 

FOR 
END 



125 ('-')',//20 nSx.13,5 x) FI 0.5,//) ) 
1600 FORM AT ( * 0 1 , 1 OX, F10 . 5 , 5X,I3) 



FAILURE RATE' , /15X, 



EXAMPLE : 



INPUT FOR THE PROGRAM 

8 . 



K= 



6 

5 

5 

5 

5 

5 

5 



.5 

5 

5 

5 

0 

0 



OUTPUT FROM THE PROGRAM 



THE INPUT IS: 
TIME= 8.0000 

K = 6 

NO. OF R.V. 



(THE N 
FAILURE 



UMBER OF DISSIMILAR LAMDA" S) 
RATE 



5 

5 

5 

5 

5 

5 



6. 5 

1. 5 

2. 5 

3. 5 
5. 0 
9.0 



FAILURE 


RATS AMOUNT 


COEFFICIENTS 




6.5 


5 








0. 15713S+G8 


0 .84332S+07 


0. 13530E+07 




0 . 2 0038E+ 06 


0 .927 01E+04 




1 . 5 


5 




0. 1 1817E + 08 




0 . 27364E+ 09 


-0 .89949Ef 03 




-0.73813E+ 06 


0 . 18586E+05 




2. 5 


5 




-0. 17601E+ 1 0 




-0. 13295E+ 1 1 


0 .823 12E+10 




0. 17341 E+ 09 


-0.1 1979E+08 




3.5 


5 




0. 25648S+ 10 




0. 18841 E+ 11 


0 .89164S+10 




0 . 25077E+ 09 


0 .334 06E+03 




5.0 


5 








-0 . 8 3480E+ 09 


-0 .52324E + 09 


-0. 16061E + 09 


o n 


-0.21713E+08 


-0 .22945E*07 




9 • r J 


-0 .39700E+ 05 


-0 . 15168E+06 


-0. 32323E+06 




-0 . 3 37 27 E+ 06 


-0 . 27860E-01 




PROBABILITY ( TIME > 


3.0)= 0.6935042 
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APPENDIX B 



This section contains a program described in Section III 
to compute the reliability of a system. 



THIS PROGRAM CALCULATES THE RELIABILITY FUNCTION: 



♦ 

♦ 

+ 

+ 

♦ 

+ 

+ 

♦ 

♦ 



PI 


F (T) 
* 


FI (T] 


( = EXP (LI ,1) 


+ EXP j 


[LI, 2) 


> . . . , 




P 2 


* 


F2 (T 


=EXP (L2 , 1 ) 


+ EXP I 


L2 , 2) 


♦ . . . « 




P 3 


* 


F3 (T 


=EXP (L3,1) 


+ EXP | 


L3 , 2) 






PI 


* 


FI (T) (=EXP (LI,1) 


♦ EXP "(LI, 2) 







ii! 



PH * FN (T) (=EXP (LN,1) + EXP (LN, 2) + +EXP (LN , MN) ) 



Ml 

> 

PROE 
TIME 
A (I, 



ARRAY FOR A 
ARRAY FOR A 
NUMBER OF E 
. TOTAL NUMBE 
A : PROBABILITY 
TIME 

K) : THE COEFFIC 
I.TH TYPE 
RATE AND K.TH CO 



LL LAMDAS. 

LL PROBABILITIES. 

XP IN EACH ROW. 

R OF ROWS. 

OF SYSTEM SURVIVAL AT TIME T. 

I ENT FOR THE CURRENT PATH, 
FAILURE 

EFFICIENT ON THE THIS POLINOMIAL 



REAL P (50) 

CALL READ 1 (IK,P ,T) 
S U M= 0 . 

DO 1 1=1. IK 

1 S UM=SU M + P (I) 

IF (ABS (SUM- 1 . ) . 
PROB A=0 . 0 

DO 2 1=1, IK 
WRITE (6 , 1009 
CALL ONECON( 

2 P ROB A= PROB A+ 
WRITE (6,601) T, 
STOP 

199 WRITE (6,198) SU 
STOP 

601 FORMAT ('I* 

*F1 0. 7) 

1009 FORMAT ('O' 

198 FORM AT ( ' 1 ' 

* EQUAL TO 
END 



GT. 1. 0 E-5) GO TO 199 



Iro.t( 

PR0*P ( 
PRO 



P(I) 

I) 



M 

,10X, 'PROBABILITY ( TIME > *,F10.5,*)= •, 



,10X, 

,10X, 

1 . 0 * 



12 , ' . ' 

•THE S 
/I OX, « 



SUBROUTINE READ 1 (IK, P 
REAL P (50) 

READ (5,499) T 



,3X, 'P(',I2, ’)=' ,F 10. 7) 

UM OF THE PROBABILITIES IS NOT 
SUM = * , FI 0. 7) 

GET INPUT 
,T) 
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WRITS (6 ,1 499) T 
READ (d, 50 0) IS 
WRITE (6,1500) IK 
READ (5,3) (P(I) ,1=1, IK) 

RETURN 

3 FORMAT (7F 1 0 . 7) 

499 FORMAT i’5X,F10. 4) 

500 FORMAT ( 5X.I5) 

14 99 FORMAT {’ 1 » ,10X, *TIME= • , FI 0. 4| 
1500 FORMAT (• 0* , 10X, *IK = *3X,I3) 
END 



SUBROUTINE ONEC ON ( PRO , T) 
REAL A (20 ,10) ,L (20) 

INTEGER NI (20) , NINIT ( 20) 
CALL READ (K,L, NINIT) 

CALL ONEDIS (K , A , NI , L) 

CALL ONEATA (K,A ,NI,L, NINIT) 
CALL CALPRO (K, A,NI,L, PRO, T 
RETURN 
END 



1 

52 

1501 



SUBROUTINE READ (K, L, N INIT ) 
REAL L ( 20) 

INTEGER NINIT (20) 

READ (5,1) K 
WRITE (6, 1501) K 
READ (5, 52) (L(I), 

RETURN 
FORM AT (5X, 15) 
FORMAT(5X,F10.5,I5) 

FORMAT ( 10 X, ' K =• ,12 ,5X, • 
* RATES. ) ' ) 

END 



NINI T (I) ,1 = 1 ,K) 



(# OF DISSIMILAR FAILURE 
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SUBROUTINE ONEDIS ( K, A , NI, L) 

REAL A (20 ,10) .L (20) 

INTEGER NI (20) 

K K=1 

A (KK , 1) =1 . 

NI (KK) = 1 
JJ=1 

IF(JJ.EQ.KK) GO TO 6 

A (KK, 1) =A (KK, 1) *L(JJ) /(L(JJ) -L (KK) ) 

IF(JJ.SQ. K) GO TO 9 

JJ=JJ+1 

GO TO 4 

IF (KK. EQ. K ) GO TO 3 

KK=KK+1 

GO TO 2 

RETURN 

END 



BEGIN TO ADD ONE AT A TIME 



SUBROUTINE ONEATAjK, A ,NI, L, NINIT) 
REAL L(20) , A (20 .10) 

INTEGER NINIT (20) , NI ( 20) 

IE= 1 

IF (NINIT (IE) . EQ.NI (IE) ) GO TO 99 
NI (IE) =NI (IE) +1 



UPDATE IE 



40 



non non 



12 



J =0 
NN= 



N=NI (IS) - J 

A (IE.NN) =L (IE) *A (I E.NN-1) /FLOAT (NN-1) 

IF (NN. EQ. 2) GO TO 20 

J=J+1 



GO TO 12 
20 CONTINUE 

DO 101 1=1, K 
IF(I.EQ.IE) GO TO 101 
A (IE, 1) =A (I 
I 

FACT 

NKK = NI (I) 

DO 102 11=2, NKK 
?ACT=FACT <11-1 
l)+A( 
CONTINUE 



24 



l (I E, 1 )=A (±E , 1 ) + A(I,1)*L 
:F(Nljl) . LT. 2) GO TO 101 
ACT — 1 « 



1) + A(I ,1) *L(I)/ (L(I) -L(IE) ) 



102 
1 01 



A (IE, 1) =A (IE, 1) + A(I ,11) *^ACT* L (IE) /(L (I) -L (IS) ) **II 

cor - 

CONTINUE 

UPDATE I.NE.IE 



1=1 

21 IF (I. EQ . IE) GO TO 26 
A (I, NX ' 



,NI (I) )=A (I, NI(I) ) *L(IE)/(L(IE)-L(I)) 
IF (Nl jl) .EQ. 1) GO TO 26 



NNI = NI (I) -J 

A(I,NNI) = (L (IE) *A(I.NNI) -FLOAT (NNI) *A (I, NNI + 1 ) ) 
* /(L(IS)-t(I) 

IF (J.EQ. (NI (I) - 1) ) GO TO 26 
J=J+ 1 
GO TO 24 

26 IF(I.EQ.K) GO TO 32 
1 = 1 + 1 
GO TO 21 
99 CONTINUE 

IF (IE . EQ. K) GO TO 35 
IS=I E+1 
GO TO 32 
35 RETURN 
END 
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SUBROUTINE CALPRO (K,A ,NI, L,PRO,T) 
REAL L (20) , A (20 , 10) 

INTEGER Nl ( 20) 

PRO= 0 . 

DO 40 1=1, X 
S U M= 0 . 

NKK=NI (I) 

DO 4 1 J=1, NKK 
S UM= SUM+ A (I, J) *T*< J-1) 

CONTINUE 
TTT=-L(I) *T 
PRO= PRO+S UM*EXP (TTT) 

CONTINUE 
WRITE (6,1010) 

DO 1 000 1=1 , K 



NNN= NI (I) 

WRITER 6^1 00 1 ) L(I) ,NI (I) , (A(I,J) ,J=1,NNN) 

10 01 FORMAT (10X,F10. 5,5X,I3,5X, ,2j5(E12.5 r 5X) ,/30X) ) 

1010 FORM AT ( 15 X, 'LAM DA ' ,5 X, 'NI* ,d X, ' COEFFICIENTS * ,/14X, 



1000 



RETURN 



135 ('-')) 

END 
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INPUT FOR THE PROGRAM 



TIME 15. 

#RO WS 5 

.2 .2 

1.ST 3 

. 1 
.8 
. 4 

2. ND 3 

. 2 
. 5 
. 6 

3. RD 3 

. 1 
. 5 
. 3 

4. TH 2 

. 1 
. 2 

5. TH 3 

. 1 
. 2 
. 4 



OUTPUT FROM THE PROGRAM : 



TIME* 15.0000 
IK = 5 

1. P (1) = 0.2 

K =3 (# OF DISSIMILAR FAILURE RATES.) 

LAMDA NI COEFFICIENTS 

0.1 0 3 0. 11294E + 01 0.SS 342S-01 0.87Q75E-02 

0.80 2 0.18742E-01 0.23324E-02 

0.40 1 -0.14815S+00 

2. P ( 2) = 0 . 2 

K =3 (# OF DISSIMILAR FAILURE RATES.) 

LAMDA NI COEFFICIENTS 

0.20 2 -0. 25077S + 02 0.23 148E + 01 

0.50 4 -0. 13017E+03 0.21 333E+02 -0.22222E-01 0.55556E-01 

0.60 1 0. 15625E + 0 3 

3. ? (3) = 0.2 

K =3 (# OF DISSIMILAR FAILURE RATES.) 

LAMDA NI COEFFICIENTS 



0.1 0 3 0. 7553 1E + 0 1 -0. 5 4932E+00 0.27466E-01 

0.50 4 -0.89945E+01 -0. 9 9536E + 00 -0.42847E-01 -0.73242S-03 

0.30 2 0.24414E+01 -0.14648E+01 

4. P( 4 ) = 0.2 

X =2 j* OF DISSIMILAR FAILURE RATES.) 

LAMDA NI COEFFICIENTS 



0.10 3 0.32000E+02 -0.16000E+01 0.40000E-01 

0.20 3 - 0. 3 10 00 E+ 02 -0.14 000E + 01 -0.20000E-01 

5. P ( 5) = 0. 2 

K =3 ( # OF DISSIMILAR FAILURE RATES.) 

LAMDA NI COEFFICIENTS 

0.10 2 -0. 15170E+ 03 0.37926E+01 

0.20 4 0. 14400E+03 0.12800E+02 0.32000E+00 0.10667E-01 

0.40 3 0.87037E+01 0.51851E+00 0.83889E-02 

PROBABILITY ( TIME > 15.0 )= 0.9111188 
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APPENDIX C 



INTRODUCTION 



This section consists of a computer program to compute 
system reliabilities as described in Section III. 

Again for simplicity, ve will go rhruogh an example. 

The structure of the system is 

spare ( Ai ) 

Ai 





• A2 

Components 1, 2 and 3 have lives exponentially distributed 

with failure rates Ai ,Xz t A 3 respectively. Also we have a 
spare for component 1 which is identical to component 1 . 

Using the shorthand approach, the system life would b<= 
determined as follows 
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Xi /<Xi+X 2 +X 3 ) EXP{X 2 +X 3 } 

EXP {\+\+X 3 } +^Xz / CVi+X 2 +X^ EXP CXi+Xa) 

Xi/iX^+Xj X 3 /^ 1 +X 2 +X 3 ) ZERO 

exp {x 1 +x 2 +x 3 } +-x 3 /{x 1 +a,+x 5 )-zero 

x 2 /V\n 3 ) Xi/tx^^j-sxp^+Xg} 

^ V 'EXP{X 1 +^ 

^ X 3 / {X+XJ-Z2R0 



Seme definitions are necessary before describing the pro- 
gra m 

NOD E: Each node represents an exponentially distrubuted 
random variable with a certain failure rate. The number for 
any node can be chosen arbitrarily but car. not be used more 
than once. 

ARC : Each arc originates at a node and leads to another 
node. Only one arc can terminate at a given node. Arc num- 
bers can be chosen arbitrarily. If n is the total number of 
nodes in a system life then there will be n- 1 arcs in this 
sys tern. 




44 



BACK POINTER LIST (IP 3) : Each node has a back pointer. A 

back pointer is an arc number which shows which arc connects 



the node to the tree. 


If the pointer is zero, then the 


related node is the root 


of the tree. 



NODE NO BACK 


POINTER IP B ( I) 


1 


0 


2 


1 


3 


2 


4 


3 


5 


4 


6 


5 


7 


6 


8 


7 


9 


8 



(Here node 1 is the root of the ores.) 

NODE CODE LIST: As we mentioned, each node represents an 
exponential lifetime, for simplicity we can use some integer 
code numbers for each failure rate. 
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Code No. 


Related Failure Rate-L(I) 


0 


(distribution ZERO) 


1 


Xi + X 2 +x 3 




^•1 + 13 


3 


Xi + Xz 



Here only the code number 0 is not arbitrary and 0 can be 
used for the ZERO distribution. The others can be picked ou- 
arbitrarily. 

ARC C ODE LIST: This list is similar to the node code list. 



Arc code numbers : 


represent probabilities. 


Code No. 


Prob ability- PA (I) 


1 


Xi /(x x +a 2 +a 3 ) 


2 


X 2 /{Xi +A 2 +A 3 ) 


3 


X 3 /{Ai*A 2 +A 3 ) 


4 


A 1 / ( A 1 * A 3 ) 


5 


A 3 / ( A 1 + A 3 ) 



ARC ORIGIN LIST : 


Each arc has an origin node and a ter mi- 



nal node 



In the program we need 



to use only the origin 



Arc No 



Arc N< 


l . Node 10(1) 


1 


1 


2 


1 


3 


1 


4 


2 


5 


2 


6 


2 


7 


3 


8 


3 


LAST POINT 


NODE LIST (LP) : This list indicates the nodes on 



the end of each path. There is no necessary order in the 



list. 

LAST POINT 


DEPTH LIST (IPD) : This list indicates the number 



of arcs from last node to the root of the tree. 



I 1 < 


, ?. Node List LP (I) L. P. Depth List IPD (I) 


1 


5 2 


2 


6 2 


3 


7 2 


4 


8 2 


5 


9 2 


6 


4 1 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



THIS PROGRAM COMPUTES THE RELIABILITY CF A SYSTEM 
WHICH HAS A SURVIVAL FUNCTION LIKE BELOW: 

? (T) = 

PI * FI (T) ( = EXP (LI , 1) +EX? (LI. 2) + + EXP(L1,M1)) 

+ P2 * ?2(T) (=EXP(L2, 1) *EXP(L2,2) + + EX?(L2,M2) 

+ P3 * FI (T) (=EXP (L3, 1) ► EXP(L3, 2) «• +EXP(L3,K3 



+ PI * FI (T) <=EXP (LI, 1) +EXP (LI, 2) + + 3XP (LI , MI) ) 

+ 

«■ 



+ PN * FN (T) (=EXP (LN f 1) +EXP (LN f 2> + +EXP (LN f MN) ) 



IPB (I) 
IPA (K) 
IO(K) 

I CL (I) 

EM! 

IPD () 

s° 

M 

NLA ST 
N DIF 
MDIF 
NPATH () 
MPATH () 
I M 
IP 

PROBA 

PRO 

P 

LCODE () 
LL() 

IK 



VARIABLES 
BACK 



Sift 

i (,) 

T 



T() 



IN PRO 
POINT E 
PROBABILITY 
ORIGIN NODE 
FAILURE RAT 
DIFFERENT P 
LAST POINT 
LAST POINT 
FAILURE RAT 
NUMBER OF N 
NUMBER OF A 
NUMBER OF L 
NUMBER OF D 
NUMBER OF D 
NODE LIST I 
ARC LIST IN 
THE NUMBER 
THE DEPTH 0 
THE CURRENT 
PROE. OF T.H 
THE CURRENT 
FAILURE RAT 
FAILURE RAT 
9 OF DIS. ? 
NUMBER DF L 
AUXILAURY A 
COEFFICIENT 
TIME 



GRAM : 

R OF THE I.TH NODE 
CODE OF THE K.TH ARC 
OF K.TH ARC 

E CODE OF THE I.TH NODE 
ROBAEILI T Y LIST 
NODE LIS I 
NODE DEPTH LIST 
E LIST 

ODES IN TREE 
PCS IN TREE 
AST POINT NODES 
IFFERENT FAILURE RATES 
EFFERENT PROBABILITIES 
N A PARTICULAR PATH 
A PARTICULAR PATH 
OF THE CURRENT LAST POINT NODE 
F THE CURRENT LAST POINT NODE 
VALUE DF THE SURVIVAL FUNCTION 
E CURRENT PATH SURVIVAL AT T 
PATH PROBABILITY 
E CODES IN CURRENT PATH 
E LIST IN CURRENT PATH 
AILURS RATES IN CURRENT PATH 
IFES RELATED TO THE FAILURE RATE 
RRAY 
S FDR 



FOR N I ( ) 
CURRENT PATH 



MAIN PROGRAM 

REAL L (20) , PA (2 0) , LL ( 1 1) , PRO, P. ?ROBA,T 
INTEGER IPB (100) ,ICL (100) , 10 (1 )0) , LP (50) , I PD (5 0) , 
♦NPATH (11) , MPATH (10) f SI (11) ,I,IX,IPA(100) 

PROB A=0 . 0 

CALL RE AD (LP , L, PA, IP3 ,ICL, 10, IPD, NLAST ,T, IPA) 
WRITE (6,1013) 

DO 1 1=1, NLAST 
I M=L P (I) 

IP=IPD (I) 

NPATH (1 ) = IK 
DO 9 J=1,I? 

NPATH (J + U = IO (IP3 (IM) ) 

MPATH (J) =I?B (IM) 

I M=N?AT H (J+1) 

I CONTINUE 
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CALL PASS (PA , L, ICL , NPAT H , M P A TH , I? , LL , HI r IK, P,IPA) 

WRITS (6, 1590) 1,1,? 

CALL 0NEC0N(PR0,T, LL,NI,IK) 

PROBA= PROS A+P* PRO 
1 CONTINUE 

WRITE (6,601) T, PRO BA 
STOP 

501 FORMAT (' 0’ , 10X, 'PROBA BILIIY ( TIME > ',710.5,')= ', 
♦F10.7) 

1013 FORMAT (10X, 'BEGIN TO CALCULATION' ,/11X,20('-') ) 

1990 FORM AT ( ' 0 ' ,10X,I2,'.' ,5X, 'P(' ,12,') = ' , F 1 0 . 7) 

END 

GET INPUT 

SUBROUTINE READ (LP,L, PA, IP B , IC L , 10 , IPD, NL A ST , T , I PA) 
REAL L ( 20 ) , PA (2 0) 

INTEGER IP3 (100) ,ICL ( 100) ,IO (190) , L? (50 ) , I ?D (5 0) , 

READ TIME AND # OF NODE 



READ (5, 1) T , N 



READ BACK POINTSRLIST, FAILURE 

RATE CODE LIST 



READ (5,3) (IPB(I) , ICL(I) ,1 = 1 ,N) 

M = N- 1 

- READ ARC ORIGIN LIST AND ARC 
PROBABILITY CODS LIST 

READ (5,4) (I O (K ) , I PA (K) ,K=1,M) 

READ # OF PATH AND # OF DIFFERENT 
FAILURE RATE AND # OF DIFFERENT 
ARC PROBABILITIES 

READ (5,5) NL AST ,MD IF ,MDIF 

READ LAST POINT LIST (L?) AND 
LAST POINT DEPTH LIST (IPD) 

READ (5,6) (LP(I) , IPD (I) ,I=1,NLAST) 

READ FAILURE RATES 

READ (5,7) (L (I) ,1= 1, NDI?) 

READ (5,7) (PA(I) ,1 = 1 ,MDI?) 

RETURN 

1 FORMAT(5X,F10.3 ,15) 

3 FORMAT (5X, 1015) 

4 FORMAT (5X, 1015) 

5 FORMAT (5X, 315) 

6 FORMAT(5X,10I5) 

7 FORMAT (5X,5F10. 5) 

END 

THIS SUBROUTINE COMPUTES NECESSARY LISTS TO 
CALCULATE THE SURVIVAL FUNCTION OF A CURRENT PATH 



SUBROUTINE PASS (P A , L , ICL, N PATH , MP ATH , IP , LL , NI , IK, ?, 
*IPA) 

REAL PA (20) ,LL (11) ,L (20) 

INTEGER NPATH (1 1) , MPATH (1 0) , L CODS (10) ,NI(10) , I CL (100) , 
♦IP A (100) 

P= 1. 

DO 1 1=1,1? 
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4 

13 



1 P=P*PA(IPA (MPATH (I) ) ) 

IPP=IP+1 

11=1 

I? (I CL (NPATH (1) ) . 2 Q. 0) 11 = 2 

DO 2 I=II,IPP 

2 LCODE(I) =ICL (NPATH (I)) 

I K=0 

DO 3 I=II ,I?P 
I F (LCODE ( I) . E3. 0) GO TO 3 
IK=IK+1 

LL (IK) =L (LCODE (I) ) 

LLL= 1 
JJ=I+1 

IF (JJ.GT. IPP) GO TO 13 
DO 4 J=JJ,IPP 
IF (LCODE (I) . NE. LCODE (J) ) 
LLL=LLL+1 
LCODE ( J)=0 
CONTINUE 



GO TO 4 



NI (IK) =LLL 
LCODE (I) =0 
CONTINUE 



RETURN 
END C 

THIS SUBROUTINE CONTROLS THE 
COMPUTING OF THE SURVIVAL FUNCTION OF THE CURRENT PATH 

SUBROUTINE ONECON(PRO ,T,L f NINII , K) 

REAL- A ( 1 1 -10) .L (11) 

INTEGER NI (11) , NINIT( 11) 

CALL ONEDIS (K,A,NI,L) 

CALL ONEAT A (K, A ,NI ,L, NINIT) 

CALL CALPRO (K,A,NI,L, PRO, T) 

RETURN 

END 



THIS ROUTINE COMPUTES 
AS EACH IS DISSIMILAR 



’HE COEFFICIENTS 



SUBROUTINE ON ED IS { K, A , NI , L) 

V" 



REAL A (11 ,10) f L (11) 

INTEGER NI (20) 

KK=1 

2 A (KK. 1) =1. 

NI (KK) = 1 
JJ=1 

4 IF (JJ.SQ. KK) GO TO 6 

A (KK f 1) =A (KK, 1) *L( JJ) /(L( JJ) - L(KK) ) 

6 IF (JJ.EQ. K) GO TO 9 
JJ=JJ + 1 
GO TO 4 

9 IF (KK. EQ. K ) GO TO 3 
KK=KK + 1 
GO TO 2 

8 RETURN 
END 



BEGIN TO ADD ONE AT A TIME 



SUBROUTINE ONEATAJK, A,NI, L, NINIT) 

, 1 * 1 ( 11 ) 



32 



REAL L ( 1 1 ) ,A (11 ,10 
INTEGER NINIT (1 1) 

I E= 1 

IF (NINIT (IE) . EQ.NI (IE) ) 
NI (IE) =NI (IE) +1 



GO TO 99 



UPDATE IE 
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12 



J = 0 

M N=NI (IS) -J 
AjfIE,NN) =L (IE) *A (I E.NN-1) /FLOAT (Nil- 1 ) 

IF (NN. EQ. 2) GO Td 26 
J=J+ 1 
GO TO 12 
20 CONTINUE 

DO 101 1=1. K 
I F (I . EQ . I sf GO TO 10 1 

AlrfE, lf=A (IE ,1 ) +A(I,1)*L(I)/ (L(I)-L(IE)) 

IFfNI (I) . LT. 2) GO TO 101 

FACT=1. 

NKK = NI (I) 

DO 102 11=2, NKK 
FACT=FACT <11-1 



1 02 
1 01 



.< FACT=FACT <11-11 

A (IE, 1) =A (IE, 1) + A(I ,11) * rACT*L (IE) / (L (I) -L (IE) ) **II 
CONTINUE 



CONTINUE 



UPDATE I.NE.IE 



1 = 1 



21 IF (I . EQ .IE) GO TO 25 

A (I, NI (I) )=A(I, NI (I) ) *L (IE) / ( L ( IE) -L (I) ) 

IF (N I jl) .EQ. 1) GO TO 26 

24 NNI = NI (I) -J 

A (I , NNI) = (L (IE) * A (I, NNI) - FLOAT (NNI) *A (I, NNI+1) ) 

IFfJ.EQ. (NI(l/-f! )* L io ) IO 26 
J=J+1 
GO TO 24 



26 


IF (I.EQ.K) 


GO TO 32 




1=1+1 
GO TO 21 


• 


99 


CONTINUE 






IF (IE . EQ. K) 
IE=IS+1 
GO TO 32 


GO TO 35 


35 


RETURN 

END 





CALCULATES THE PROBABILITY FOR THE CURRENT PATH 

SUBROUTINE C ALP RO ( K, A , NI, L,PRO, T) 

- ... A / V 

(IT) 



REAL L (11 ) , A (11 , 10) 
INTEGER NI ' 
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PRO=0, 

DO 40 1=1, K 
SUM=0. 

NKK=NI (I) 

DO 4 1 J=1, NKK 
S UH=SUM+ A (I, J) *T*< J-1) 
CONTINUE 
TTT=-L (I) *T 
IO+S 



FRO= PRO+S UM*EXP (TTT) 

40 CONTINUE 

WRITE (6,1011) 

DO 1000 1=1, K 
NNN= NI (I) 

WRITE (6,1010) L (I) , NI (I) , ( A (I , J) , J= 1 , NNN) 

1000 CONTINUE 
RETURN 

10 10 FORMAT (’O’ , 10X, F10.7, 5X, 1 2 , 5X , 2 (5 ( El 2. 5 , 2X) ./2 5X) ) 
10 1 1 FORMAT (15X, 'LAM DA' ,7X NI' ,5X, ' COEFFICIENTS ' ,/10X, 
135 (•-•) ) 

END 
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INPUT DATA FOR THE PROGRAM: 



T.N 

N<5. 1 

NO. 6 
NO. 1 1 
AH. 1 
AR. 6 
AR. 1 1 
NL, MN 
LPI PD 
LP6PD 
LAM DA 
PA' S 



2 . 

0 

5 

10 

1 

3 

6 
6 
3 

13 

2 . 0 

. 25 



1 

1 

3 

1 

2 

2 

3 

3 

3 



I 3 
1 

6 

I I 
1 
4 
7 
3 
9 



1 

2 

2 

2 

3 

3 



2 

7 

12 

2 

5 



10 



1 

3 

2 

1 

1 



3 

8 

2 

5 



1 1 



3 

3 

2 

2 



4 

9 

3 

6 



12 



5 

75 



1.5 

1.0 



OUTPUT FROM THE PROGRAM 



BEGIN TO CALCULATION 

1 . ? ( 1)= 0. 0625000 

LAM DA NI COEFFICIENTS 

1.5 2 - 0. 80 000E + 0 2 0. 24000E+02 

2.0 2 0.81000E+02 0. 18000E+02 

2. 'P( 2) = 0.0 46 875 0 

LAM DA NI COEFFICIENTS 



1.5 2 -0. 51 200E+03 0. 96000E+02 

2.0 3 0. 51300E+03 0.16200E+03 

3 . ? ( 3)= 0. 1406250 

LAM DA NI COEFFICIENTS 



0.5 2 -0. 63578E-06 0. 11852E + 01 

2.0 3 0. 10000E+0 1 0. 81481 E+00 

4. p( 4)= 0.0 46 8750 

LAM DA NI COEFFICIENTS 



1.5 2 -0. 51 200E+0 3 0. 96000E+02 

2.0 3 0. 51 300E + 0 3 0. 16200E+03 

5. P ( 5 ) = 0. 1406250 

LAM DA NI COEFFICIENTS 



0.5 2 -0. 63578E-06 0. 1 1852E+01 

2.0 3 0. 10 000E + 0 1 0. 81481 E+00 

6. ? ( 6) = 0.5625000 

lam da ni Coefficients 



0. 1 8000E+-0 2 



0 . 22222E+-00 



0 . 18000E+02 



0. 22222E+00 



0.5 2 0. 59259S+ CO 0. 8 8889E + 00 

2.0 2 0.40741E+00 0. 22222E + 00 

PROBABILITY ( TIME > 2.0 )= 0.0173515 



1 

2 

1 

1 



3 
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APPENDIX D 



This section gives t he program for the simulation of a 
system's life mentioned in Section III. The program uses a 
crude Monte Carlo simulation procedure. 



THIS PROGRAM SIMULATES A SYSTEM HAVING 
THE RELIABILITY FUNCTION: 



+ 

+ 

+ 

+ 

+ 

+ 

♦ 



P (T) = 



PI 


* 


FI 1 


[T) 1 


r=EXP' 


(LI, 


1 ] 


+ EXP 


(LI , 2) 


+ . . . , 




(LI ,M1) 


P2 


* 


F 2 1 


T) 1 


= EXP i 


L2, 


1 ' 


+ EXP 


(L2,2 


+ . . . , 




L2,M2 


P3 


* 


F3 | 


T { 


= EX P ■ 


L3, 


1 ( 


♦ EXP 


(L3,2 


+ . . . , 




L3,M3 



PI * FI (T) ( = EXP (LI, 1 ) + EXP (LI, 2) + +EXP (LI, MI) ) 



+ 

+ 



PN * 




PROB 

TIME 



FN (T) (=EX P (LN, 1 ) + EXP (LN , 2 ) ♦ +EX? (LN,MN) ) 

: ARRAY FOR ALL LAMDAS. 

: ARRAY FOR ALL PROBABILITIES. 

: NUMBER OF EX? IN EACH RON. 

: FIST NUMBER OF LAMDAS IN THIS ROW. 

: TOTAL NUMBER OF ROWS. 

; PROBABILITY OF SYSTEM SURVIVAL AT TIME T. 

: TIME 



REAL L (50 0) ,P(S0) ,?ROB 

INTEGER MI (50) , MT(50) , MIDO , MIUP , I r J, N 

IX=4 56378 

CALL READ (L,N, TIME, MI ,ST, P) 

CHECK FOR SUM OF P'S EQUAL TO 1.0 



TOT=0. 

DO 150 1=1, N 
150 TOT=TOT+P (I) 

IF (A3S (TOT-1. ). GT. 1. E-5) STOP 

CALL CONTRO (L,N ,TIME, MI ,MT , P, IX , P ROB) 

CALL OUTPUT (L,N ,TIME, MI,MT,P, FROB) 

STOP 

END 



SUBROUTINE READ 



C 



SUBROUTINE READ 
REAL L (500) ,P (5 
INTEGER MI (50), 



MT (53) 



IME, M I, MT, 



P) 
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GET THE TIME 

READ (5,505) TIME 
GET THE NUMBER OF ROWS 
READ (5,501) N 
GET ALL P'S 

READ (5,504) (P(I),I=1,N) 

GET ALL LAMDAS IN THE ORDER OF ROW BY ROW 

MIDO=1 
MI UP =0 

DO 100 J= 1 , N 
IF(J.NE.I) MI DO =MI UP+ 1 
MT (J) =MIDO 
READ (5, 520) MI(J) 

MIUP= MIUP+MI (J) 

READ (5, 503) (L ( I) , I = M IDO, MIUP) 

1 00 CONTINUE 
RETURN 

501 FORMAT (51,15) 

503 FORMAT (5F10. 5) 

504 FORMAT (5F10. 7 

505 FORMAT(5X,F10.3) 

520 FORMAT (5X, 15) 

END 



SUBROUTINE 

SUBROUTINE CONTRO (L, N ,TIME , MI , MT , P , IX , PRO) 
REAL L ( 50 0) - P (5 01 
INTEGER MI (50) , MT(50) 

SUM=0. 

DO 1 1=1 -N 
NN = 1 000000*? (I) 

CALL SIMULA (NN,L,M 1(1) ,MT(I) , X, IX, TIME) 

1 SUM=SUM+X 
PRO= SUM/1 C00000 . 

RETURN 

END 

SUBROUTINE FOR SIMULATION 

SUBROUTINE SIMULA(NN, L, M, MTT, X, IX, TIME) 

REAL L ( 50 0) ,RN(50) 

INTEGER M,I,J,IX 
X = 0 

MMT= MTT+M-1 

DO 1 1 1=1, NN 
TEST=0 . 

CALL LEXPN (IX, RN,M ,16307,0) 

JJ=0 

DO 111 J=MTT , MMT 
JJ=JJ+ 1 

111 TEST=TEST+RN (JJ)/L(J) 

11 IF (TEST. GE. TIME) X = 1.0 + Z 
RETURN 
END 

SUBROUTINE 



SUBROUTINE OUTPUT (L, N , TIM E, MI , M T , ? , PRO) 
REAL L (500) ,? (50) 

INTEGER MI (50) , MT(50) 



CONTRO 



OUTPUT 
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PRINT OUT ALL THE SYSTEM 



WRITE (6,600) 
MIUP=MI (1) 



HHI Sg (L(I) 



,1= 1 , M IU P) 



MIDO=MT (I) 

MIUP=MT(l[ + MI(I) -1 
WRITE (6,60 2) P (I) , (L (J) , J=MID0 , H IUP) 
102 CONTINUE 

WRITE (6 ,603) TIME, PRO 
RETURN 



FORMAT STATEMENTS 



FORMAT (' 1 1 , 5X, 1 THE SYSTEM IS :•) 
FORMAT ( , 0 , ,5X, , F(T)=' ,F10 . 7 , ' * 

I (F10.4,3X) ,5 (//20X,5fF10. 4,3X) ) ) 
FORMAT ( '0 ' , 9X, ' + ',FT0.7,« * *,»EXE 



6 CO FORMATS 1 ' , 5X,' THE SYSTEM IS 
601 FORMAT( '0 * ,5X,' F ‘ " 

1 " “““ 

6 C 2 rwnoaii '« ' i 7Aj ■ t r r i u • / • ' 

10 . 4 , 3 X) , 5 (//20X,5(F10 .4,3X) ) ) 

603 FORMAT* '0 • f 5X,' PROBABILITY OF SYSTEM SURVIVAL AT 
*T=',F10.5, f IS 1 , F 10 . 7) 

END 



EXP L = ' , 5 
EXP L = ' ,5 (FI 



INPUT FOR THE PROGRAM 



TIME 
#RO W 

. 2 

1. ST 

. 1 
. 8 

2. ND 

. 2 
. 5 

3. PD 

. 1 
.5 

4. TH 

. 1 
. 2 

5. TH 

. 1 
. 2 



15 

5 

6 
7 
9 
6 
9 



.2 

.8 

.2 



.5 

.5 



.2 

.2 



.2 

.4 

.5 



.3 

.5 

.1 



.4 

.4 



2 

1 



.2 
. 1 

. 5 

. 1 

. 2 

.2 



OUTPUT FROM THE PROGRAM 



THE SYSTEM IS 
F (T) = 0.2 * EXP 



+ 

+ 

+ 



0, 

0. 

0. 

0. 



2 
2 
2 
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SURVIVAL AT 


T = 


15.00000 


' IS 



0.911241 
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APPENDIX E 



This section reviews some notions which are found in the 
references for this paper. 

E. 1 Redundant systems with exponentially lived compo- 
nents : 



In reliability analysis, 
describe a set of components o 
sion. A system is redundant if, 
components involved act as back 

A rough definition might be 
tant if the failure of any one 
failure of the system, and is 
its components can fail without 
Thus redundant systems have the 
fer damage through the failure 
and still survive. (SSARY'Ref.1 

E. 2 First failure in a set 
ner. ts : 



he term system is used to 
ganized to perform some mis- 
in some fashion, some of the 
up for other components. 



that a 


system is 


not redun- 


of its 


components 


causes the 


redunda 


nt if one 


or more of 



causing the system to fail, 
property that they car. suf- 
of some of their components 

]) 

of exponentially lived ccmpo- 
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We have n components, each indeoende.nt from the others, 
we want to compute the probability that the j.zh component 
fails first. 



P(j.th component fails first) =P(T <T. ,\/i *i*j) 

J i 



-A .s 



cr 



P ( 1 . >T- , v i,U#j|I =s) i e J ds 



o' 



[n F. (s| ) ]f (s| ds 

i^j 1 J 

-\, s "V 
[ n e 1 ]A.e J ds 

its J 

n n 



3 A/< I X. ) [ l ] a 3s 

J i=l ’ °l i=l 1 
n 

= X/( l X. ). 

J i 1 

E. 3 Degeneracy at zero (Zero Distribution): 

Let ZERO be the name for the distribution of a random 
variable that is degenerate at zero. 

If P(T o =0)=1, then we say rhat has the distribution 

Z ER 0 , or 



T ^ ZERO . 
o 
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